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We study the crossover from the one-dimensional to the two-dimensional Hubbard model in the 
photoemission spectra of weakly coupled chains. The chains with on-site repulsion are treated using 
the spin-charge factorized wave function, that is known to provide an essentially exact description 
of the chain in the strong coupling limit. The hoppings between the chains are considered as a 
perturbation. We calculate the dynamical spectral function at all energies in the random-phase 
approximation, by resuming an infinite set of diagrams. Even though the hoppings drive the system 
from a fractionalized Luttinger-liquid-like system to a Fermi-liquid-like system at low energies, 
significant characteristics of the one-dimensional system remain in the two-dimensional system. 
Furthermore, we find that introducing (frustrating) hoppings beyond the nearest neighbor one, the 
interference effects increase the energy and momentum range of the one-dimensional character. 



I. INTRODUCTION 

The Hubbard model is believed to contain most of the fundamental physics of a great variety of materials ranging 
from weak interacting metals to Mott insulators. It may also capture some of the phenomena responsible for the 
high Tc superconductivity of cuprates. In one spatial dimension the Hubbard chain is exactly solvable by the Bethe 
ansatz [1 , furthermore, the low energy properties are understood in details within the Luttinger Liquid (LL) theory 
pj. In particular, the elementary excitations turn out to be fractionalized, carrying either charge or spin quantum 
numbers, a property that is a rather generic feature in ID strongly interacting electron systems Thanks to 
the advent of bosonization the low energy physics of ID interacting electrons is now well understood including the 
computation of many observables However, despite the enormous success of these complementary approaches 

the computation of observables for arbitrary energies was obtained only in some restricted limits ^8HlO| or using some 
additional approximations |llH15j . 

In higher dimensions the physical picture is much less clear. It has fueled controversy, mainly motivated by the 
hight Tc phenomena in the layered cuprate oxides. Generically, in dimensions higher than one, excitations are not 
fractionalized, the most well-known example being the quasi-particles in a Fermi liquid (FL). A few remarkable 
experimental and model examples exist, however, where the electrons fractionalize. The most spectacular example is 
the fractional quantum Hall effect, where fractionalization of quasiparticles has been predicted theoretically |T5] and 
consequently found experimentally |17| . Fractionalization has also been theoretically proposed in electron systems 
with frustrated nearest neighbor interactions |18j . Further examples include quantum spin- liquids, where the presence 
of frustration may lead to deconfinement of the spinons in the two-dimensional system [T^. For example, the quasi-two 
dimensional triangular spin system Cs2CuCl4[2D^ has an excitation spectrum that can be described, similarly to the 
one dimensional case [21j, by a continuum originated from fractionalized pairs of spin 1/2 spinons|22]. This property 
has been verified experimentally for several quasi-one dimensional spin 1/2 systems like CPC [53], KCUF3 |24l [25] 
and copper benzoate [26'. 

Experimentally, the single-particle properties of the material are most directly measured by photoemission. The 
intensity of the extraction of the electron by photon at given energy and momentum transfer is directly proportional 
to the spectral function - the imaginary part of the one-particle Green's function. If most of the spectral weight is 
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carried by well defined spectral lines one expects excitations to be sharp, electron-like coherent modes (quasiparticles) . 
On the contrary, broad continua signal fractionalization of the electronic degrees of freedom. 

Low energy descriptions for the two dimensional case have been proposed that predict a fractionalization description 
of the low energy physics. Experimentally such low energy features are difficult to observe in the photo-emission data 
since they are obscured by resolution and noise. Therefore, it is useful to have a prediction over the full energy range 
to compare with experiments. The motivation for this work is twofold: on the one hand to provide an approximate 
spectral function valid to arbitrary energy and on the other to clarify the role of frustration in the underlying excitation. 

In this work we address the dimensional crossover from one to two dimensions in a strongly correlated electron 
system by coupling Hubbard chains within the random phase approximation (RPA). This approximation leads to a 
description of the 2D quantities spectral function in terms of the ID Greens's function of the chain. Though coupled 
one-dimensional chain tend to order at low temperature, in this work we will assume that we are at sufficiently high 
energies uj and temperatures T , typically larger than some crossover temperature Tid after which LL perturbation 
should be valid. The general expression obtained by the RPA for the two dimensional spectral function is valid for 
any values of U and filling factor, as well as for any kind of small inter-chain hopping. In particular it is possible to 
study the role of hopping in different geometries, highly frustrated cases as well as non-frustrated ones. Due to the 
lack of theoretical expressions for the spectral function in one dimensional for generic U we concentrate our study on 
the [/ — > oo limit using the results derived in [27J . The exact results obtained for the spectral function of the Hubbard 
chain in the [/ — >■ oo limit, can be extended to finite but large U and are used in the RPA to obtain the same function 
in higher dimensions. 

Previous works have dealt with the issue of coupled LL or coupled Hubbard chains. Contrasting with the LL-like 
features of decoupled chains, FL behavior is generically expected for weakly interacting systems and large inter-chain 
hopping terms. The interpolation between LL and FL regimes as the inter-chain hopping increases, as well as the 
energy scales for which each description is valid have been largely discussed. 

Using perturbative renormalization group (RG) and an RPA-like expression for the two-dimensional Green's func- 
tion, it was shown [28J, starting from a LL, that the hopping is relevant \i 6 < 1 and irrelevant if > 1, where 9 is 
the LL exponent characterizing the low frequency behavior of the density of states N (ui) ~ , (note that 9 = 
corresponds to the non-interacting case). In the first case the two-dimensional Green's function develops well define 
poles near the Fermi energy with a nonzero quasiparticle residue (Z) for non-vanishing inter-chain hoppings and in 
the second case Z vanishes. In the same direction it was pointed out that using a d = 1 -I- e expansion that the only 
weak-coupling fixed point for e > is the FL one [29] . Using a path integral formulation [3D] (like RPA) the results of 
[5S] where rederived, but it was pointed out that higher order processes could extend the FL behavior beyond 9 = 1. 
Subsequent works, using exact resummation of some infinite class of diagrams |31H33| . also corroborate this result. 
It was also shown, using bosonization, that even if long range 3D Coulomb interactions were considered the ID LL 
regime leads to a FL, for any hopping, but anomalous scaling was found in the FL phase for small hoppings |341I35| . 

The picture that FL behavior is obtained as soon as inter-chain hopping is introduced has, however, to be interpreted 
as being valid only above some finite energy scale. The introduction of inter-chain hoppings will in general lead to 
instabilities towards some possible ordered phases. The phase diagram of a system of coupled chains, including 
ordered phases, was studied in Refs. |36fHn] . e.g. The FL behavior appears for energy scales higher then the critical 
temperatures of such ordered phases. 

Moreover, for energies higher than some characteristic energy of the order of the inter-chain hopping amplitude 
(possibly renormalized by the interactions), one expects to recover LL features. Thus only for intermediate energies 
is the FL picture expected to hold. Indeed in Ref. [41 it was argued that even though the transverse hopping is 
a relevant perturbation, in the RG sense, incoherent single particle hopping between chains can lead to a LL-like 
behavior. This was confirmed in |42l H5] using exact diagonalizations (ED) and quantum Monte Carlo (QMC) since 
the incoherent part of the spectral function (SF) is less affected by inter-chain hopping, and the Drude weight is 
small compared to the incoherent weight, even for small 9. For larger 9 the hopping between chains becomes fully 
incoherent. Furthermore, considering a higher dimensional mesh of coupled LL it was shown that there are mixed 
characteristics of LL and FL jUj . 

A rather unifying picture was obtained using chain dynamical mean-field theory (CDMFT)|45|. These studies 
observe a crossover from a LL at high temperatures to a FL at low T with the coexistence of a Drude feature with 
small spectral weight and a large incoherent weight. 

Several studies also treated the case of coupled Mott insulators. The RPA approximation was used in t46j and it 
was found that for high enough hopping and small enough Coulomb coupling the Mott-Hubbard gap closes and small 
Fermi pockets appear in the Fermi surface with a finite Z. However, it was shown using CDMFT that when the gap 
closes there is a continuous FS and no pockets [45_. These results were also confirmed in |47] but it was found that 
between the Mott phase and the FS phase there is an intermediate phase where there are pockets (arcs because of 
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spectral weight inhomogeneities) . Defining the FS both by the poles and zeros of the ReG(cL', k) it was shown that 
the Luttinger Theorem [48_ is satisfied. 



The paper is structured as follows: In section [IT] we discuss the model and method used, briefly reviewing the 
RPA approach. In section |III| we present results for the spectral function at low energies where a Luttinger-liquid- 
like universal description holds and compare the results with other methods previously obtained. In section IV we 



consider the regimes of finite energies and consider finite but large U values, the infinite U limit where the spins 
are dispersionless and the half-filing Mott-insulator case. In section |V] we study the role of frustration comparing a 
square, a triangular and a fully frustrated lattices. We present some conclusions in section |VI[ Also, in Appendix 
[A| we review the method of Ref. j21] developed for the spin structure factor of the Heisenberg antiferromagnet in a 
triangular lattice and present its generalization to the electron spectral function. In Appendix [B] we briefly review 
the method used for the calculation of the spectral function for the Hubbard chain. In Appendix [C] we review the 
derivation of the RPA formulation and derive the expansion for its leading correction. This involves the knowledge of 
higher correlation functions for the Hubbard chain, which are not available at this time. 



II. MODEL AND METHOD 



This section presents the method used to obtain the spectral function of the weakly coupled Hubbard chains in 
terms of the one dimensional spectral function. In order to set the notation we write the Hamiltonian for the 2D 
Hubbard model as sum of an intra and an inter-chain term. 



U,y 



here 



x.a 



is the intra-chain contribution to the Hamiltonian of a chain parallel to the x direction and the subscript y labels the 
direction perpendicular to the chains. The hopping amplitude between the sites in the chain is denoted by t, while U 
is the usual on-site repulsion that penalizes doubly occupancy of a given site. The transverse term is given by 
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Figure 1: Direction of the hopping terms in the 2D Hubbard model. 
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where <^ labels the diflFerent inter-chain hoppings along the directions 6i displayed in Figjl] As shown, setting 
^2 — ^3 — ^ corresponds to an anisotropic square lattice, ^2 = ^'i^d = to an anisotropic triangular lattice, and 
t'^ ^ and t'2 — tg to the square lattice with diagonal hoppings. 



A. RPA and spectral function 

As briefly reviewed in Appendix [C] the single electron Green's function in the so called random phase approximation 
(RPA) is given by 

G(w,k) = [G-^{Lo,K)~t' {k)]-' (1) 

where 

i'(k) = -2^<^cos(k.^,) (2) 

i 

= —2 [t[ COS (ky) + t'2 cos [ky + kx) + t'^ cos {ky — kx)] 

is the Fourier transform of the hopping matrix. Here Gid is the Green's function of the one-dimensional system, 
assumed to be known. In this work it is calculated exactly. The Fermi momentum kj? and the QP weight are obtain 
from Eq. ([T]) requiring 

G-i(cj = 0,ki.) = (3) 
=9„G-i(w = 0,kj.). (4) 

In several works pioneered by Wen |[55^ this expression has been used to study weakly coupled Luttinger Liquids 
|30, 49_. Note that Eq. ([ij is exact for non-interacting electrons (0 = 0). 

Using Eq. ([T]) and the asymptotic form of the retarded Green's function, in the low energy limit given by bosonization 
and parameterized by 

^^K^ + ^-A (5) 



4 V is:, 

where Kc is the Luttinger parameter, it was shown [U [251 HDJ HH] that for < \ there is a nonvanishing QP weight 





i'(k) 




A 



(6) 

for an arbitrary t' ^ 0, where A is an energy cutoff and 7 is a exponent that can be explicitly computed (see chap. 
19 of |4j). Note that for the non-interacting case 6* = 0, the low energy regime of the infinite U limit of the Hubbard 
model is recovered setting 9 — 1/^ and higher values of 6 correspond to models with long range interaction. Besides 
the region 6 > 1. where no coherent mode is found at the RPA level in \A, the authors considered the regimes 
< 1/2 and 9 > \/2 for which the exponent 7 in ^ changes from positive to negative. They concluded that the 
value of the QP residue will be larger in the second region. We will see further that there is a clear physical signature 
separating these two regimes. 

As stated in the introduction t' is a relevant perturbation in the RG sense, and thus the above treatment is valid 
only for energies T,uj > Tc , where is the highest critical temperature of all the possible order phases towards which 
the system is unstable at low energy. Another energy scale is defined by Tix > Tc which separates a low energy regime 
where the pole of the Green's functions is physically perceptible [^] from an higher energy regime for which fully 
coherent 2D hopping is suppressed. For the non-interacting case (6 = 0) Tix is of the order of the interchain coupling 

t'. It has been shown [30 that for the interacting case this energy scale is reduced yielding Ti^ ^ A( 

for 9 < 1. For 9 > I this treatment leads to a vanishing Ti^', however, as also noticed in [301, higher order terms 
that consider two-particle processes define another energy scale that will overtake Ti^ and further extend beyond 
9 — 1 the region where Z ^ 0. Note that these works are only valid for arbitrarily small energies since the one 
dimensional quantities are given by bosonization and thus no predictions can be obtained for the moderate and high 
energy regimes. One of the aspects of the present work is precisely to be able to access these regions. 

Another feature of the RPA expression is that it leads to an anisotropic QP weight along the FS which vanishes for 
t'(k) = 0. This could suggest the existence of hot-spots in the FS where the ID character would be strongly manifested. 
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Figure 2: Left Panel: Values of the charge exponents /3^and fi^ as a function of the LL parameter 6. Central and Right 
Panels: Spectral function obtained using the RPA expression for t' = 0.2 and t' = —0.2 respectively and for 9 = 1/8, Vs = 1 
and Vc — 2.718. The Gray lines signal the boundaries of spin and charge continuum of the ID spectral function: u — ±Usk 
and UJ = ±Uck. The bound state (Green line) is obtained solving Eq. ([3| . Red and Blue dots displayed at tj = correspond 
to values of kx for which signRe [G {lj = 0, k^)] is respectively positive or negative. 



However, subsequent works, using exact resuniniation of some infinite class of diagrams pTI - 155] and higher dimensional 
bosonization f3S] pointed out that the vanishing Z was an artifact of the RPA and that the inclusion of higher order 
terms leads to a smoothly varying QP along the FS; this fact was also verified by DMFT calculations [50-53_. All 
these works predict a finite QP pole leading to FL like behavior for non-zero t' for the Hubbard model. 

However, the RPA expression gives a qualitative description of the crossover from 1 to 2D. In practice the use of RPA- 
like expressions has gathered a great success describing antiferromagnetic spin chains [22 , 54_ whith a good quantitative 
agreement with experiments. In electronic systems the DMFT approach, based in a large D± (dimensionality of the 
transverse dimension) expansion, obtained a good agreement for the frequency dependent interchain conductivity jSlj . 

In this work we use two approaches to compute the one dimensional spectral function. The first is only valid for 
low energies and is equivalent to the use of the asymptotic Green's function given by bosonization. It was used to 
verify the predictions referred in the last section and to understand the low energy limit of the second approach valid 
for arbitrary energies. Due to its generality it permits to vary independently the interaction strength (changing 6) as 
well as the spin and charge velocities. The second approach relies on the exact solution of the large U limit of the 
Hubbard model. This limit permits considerable simplifications and in particular a closed form for the ID spectral 
function. A detailed description of both methods in given in the following sections. 

The lowest lying excitations contributing to the spectral function of the ID Hubbard model correspond to the 
creation of a holon and a spinon (charge and spin excitations). These two quasiparticles propagate with different 
velocities and in terms of the original electrons are very complex. Even though they have a fractionalized existence 
inside the ID many-body system, when an electron is, for instance, removed from the chain (photoemission) they 
recombine. If the chains are weakly coupled one expects that the excitations travel along the transverse direction 
as "electrons". The holon and the spinon are expected to propagate coherently from one chain to the next. This 
idea was proposed in Ref. |22| in the context of an antiferromagnet in a triangular lattice. In the ID Heisenberg 
antiferromagnet the low lying excitations are two spinous. In the weak coupling regime they are assumed to propagate 
coherently from chain to chain (like a AS — 1 excitation - a magnon). It is therefore interesting to generalize the 
procedure developed in Ref. [22J, for the spin structure factor of the antiferromagnetic Heisenberg model, to the 
present case of the spectral function of the Hubbard model. This is carried out in Appendix |^ There are however 
difficulties associated with instabilities of the system resulting from the approximation used. The expression obtained 
for the spectral function is formally very similar to the one obtained within RPA (Appendix |C| as noted in ref. [22] 
for the antiferromagnet. A basic difference is that in the RPA the spectral function is defined as a complete function 
(for positive and negative energies) while in the restricted Hilbert space considered in Appendix [A] the positive and 
negative energies are associated with two functions defined separately. Due to the appearance of bound states one 
is led to a situation where the excited states have negative energies, which implies an instability of the groundstate. 
Therefore we will use in the following the RPA expression ([ij to obtain the 2D Green's function. In this context the 
bound states are interpreted as coherent modes resulting from spectral weight transfer among different energies, as 
discussed next. 
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To fix the notation we define the spectral function as 



1_ 



Sp {liJ, k) = IniG (cj, k) . 

TT 

In the Hterature it is usual to write this quantity as a sum Sp (w, k) = A{io, k) + B{lu, k) where 

A{u;, k)=J2\{f,N+l\ cl^ |0, N)\'siu; ~ Ef+' + 
f,<y- 



(7) 



(8) 



is the measured amplitude of angular resolved inverse photoemission experiments, here given in the Lehmann repre- 
sentation, and 



Biuj, k) = ^ I (/, iV - 1| cfc,. |0, N)f Sicu - < - Ef-') 



(9) 



the measured angular resolved photoemission amplitude. is the number of electrons, and/ denote the ground 
and final states respectively, the chemical potential is taken such that the ground state corresponds to zero energy so 
A{uj < 0, k) = and B{uj > 0, k) = 0. 



III. SPECTRAL FUNCTION AT LOW ENERGIES: LUTTINGER-LIQUID-LIKE REGIME 

In this section we concentrate on the low energy region that is characterized by linearized dispersions and power-law 
behavior, and study how the 2D spectral properties for low energies emerge as a function of t' and 6 within the RPA 
0. We recover some results by other authors, reviewed in the last section, and find some new features characterizing 
the different regimes. 

For low energies, and near the Fermi momentum, the spectral function of one dimensional gapless electronic systems 
can be written as a convolution of the spin and charge parts 

Sp (c., fc) (X <o<'.r " - ^,3) ^KK^.+K,,,, +5{lo + + Ql,^^,) 6k,-K.„-K,^^] (10) 

where Ki,j — 27r(i — j)/L are the momenta of the excitations, fif^- = 2TTVa{i + j)/L are the corresponding energies 
(with a = c, s and Vs and Vc are the spin and charge velocities) and their weights are explicitly given by 

a ^ F (^ + /3+ + 1) F {j + /3- + 1) 

z!F(/3+ + l) j!F(/3^ + l) • ^ ' 

The exponents /3+,/3~ and I3'^,fi^ characterize the divergence of the spectral function at the edges of the (right, -|- 
and left, — ) charge and spin continua at either the right or left Fermi points. For a Luttinger liquid with SU(2) spin 
rotation symmetry both = {^5,—!} are fixed. The charge exponents are given by 



{ff,« = {^i.^I} 



(12) 



where 9 is related with the Luttinger parameter Kc by Eq. ([5| (see also Figj2|. As we already mentioned, = for 
the noninteracting fermions, and — >■ 1/8 for J7 — > -l-oo Hubbard model. 

The particular form of the spectral function given by Eq. ( To]) was obtained in Ref. [S] in the context of the large 
U approximation of the Hubbard model. However, the described low energy structure is much more general and can 
be traced back to the conformal invariance of the (1+1)D model [55J. In the thermodynamic limit one obtains the 
well known asymptotic form, say for the right moving electrons, of the real time Green's function 



Gr(x,t> 0) ~ 

1 

{x - vsty+(^^ (x + vsty+t^i 
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(9 = (9 = 1/8 ^ = 3/4 ^ = 2 




Figure 3: Spectral function obtained by RPA formula ([T]) for different values of the LL parameter 6 and inter-chain coupling t' . 
Axes and labels are the same as for Figs. [2]- Central and Right Panels. For = we set Vc ~ «s in order to obtain the exact 
free particle result; for all other values oi 9 Vs = 1 and Vc — 2.718. The Gray lines signal the boundaries of spin and charge 
continuum. Red and Blue dots displayed at oj = correspond to values of for which signRe [G (oj = 0,kx)] is respectively 
positive or negative. For the cases where Z ^ this criterion corresponds to kx being inside or outside the Fermi surface. The 
Green line corresponds to the bound states obtained solving Eq. ([3|. 



that can be directly obtained by bosonization techniques. 

With the ID Green's function computed with (10 1 we used Eq. ([ij to obtain the 2D spectral function S'p(w,k) — 
— ^ImG'(w,k) for a fixed value of <'(k) = t' . In Figure |2]- (Central and Right panels) we show the typical results 
obtained here. The bound states were found solving Eq. (|3| outside the spin and charge continua. The FS was 
determined for the values of k for which ReG'(a;,k) changes sign. Figure |3] displays the main results of this section. 
The spectral function is shown for different values of the LL parameter 6 and inter-chain coupling t' . For = we set 
Vc ~ Vs in order to obtain the exact free particle result; for all other values of 9, fixed values of the spin and charge 
velocities were used for the physical case Vs < Vc- Fig. [i] shows the values of the QP residue as a function of t' for 
different values of 0. The error bars are due to the discreteness of the k values: for each value of t', fc+ and k^ were 
determined on each side of the FS. For these values the bound state equation was solved in order to find ~ 0; 
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Figure 4: Quasiparticle residue as a function of t' for different values of 6: = 1/16 (Orange), 6 = 1/8 (Blue), 6 = 1/2 (Gray) 
and 6 = 1 (Black). For each value of t' , and fcj were determined on each side of the FS. For these values the bound state 
equation was solved in order to find 0J+/- — 0, ^+/- was computed using Eq. Q. 

was then computed using Eq. 
As a general feature, we note the change from incoherent regions, arising from the spin-charge separation in the 
ID case {t' — 0), to the sharply defined coherent excitations as t' increases. For the ID case the spectral function 
is strictly zero outside the ID continuum, delimited by the spin and charge velocities (see |3] upper- left panel). The 
interchain coupling t' favors the appearance of sharp coherent features not only outside the ID continuum, where 
they correspond to poles of the 2D Green's function, but also within the ID continuum where the spectral weight also 
tends to concentrate. For 6 < 1 and small positive i'(~ 0.1) there is a considerable transfer of spectral weight to the 
spin (charge) branches for a; < (a; > 0) . For negative t' the spin and charge roles are interchanged (see Fig. 2]- 
Central and Right Panels). The critical value of 6* = 1, predicted by several authors 11112111301111], is found such that 
for < 1 a bound state appears for t' crossing a; = at kp{t') ^ kp{t' — 0) changing the position of the Fermi 
surface and resulting in a non- vanishing QP weight Z. For 6 > I Fig. [s] shows that for small values of t' the bound 
state still appears. However, since it does not cross w = 0, it is unable to drive the system to a FL like behavior. 
After some critical t' is reached the bound state crosses twice the lu = line creating a hole pocket. Note, however, 
that in this regime the RPA is probably out of its domain of validity and this last feature is probably an artifact. 
In figure Fig. [sjwe show the evolution of the quasiparticle residue as a function of theta. For the region < 1/2 a 
damped mode is observed when the coherent mode enters the charge continuum. For 6 > 1/2 this feature disappears 
and the coherent mode is defiected to a; = and loses all its spectral weight before entering in the continuum. This 
feature clearly differentiates both regimes. The QP residue as a function of t' is shown in Fig. |4j The large error bars 
obtained due to the discreteness of the values of k prevent a clear fit. 



In this section we use the spectral function obtained for large U [27], together with the RPA expression ([T|), to 
compute the finite energy spectral function, for weakly coupled Hubbard chains. The results presented here generalize 
to finite energies the ones obtained in the previous section for systems that can be well described by an Hubbard like 
Hamiltonian, with relatively large onsite repulsion (U/t > 6). 

It has been shown that in the U ^ oo limit the eigenstates of the Hubbard chain can be written as a product of 
a spinless free fermion and a squeezed spin wave functions |56l [57] . In subsequent works |S| [22 this factorized form 
was used to write the spectral function as a convolution over the spin and the fermionic parts (see Appendix IbI . The 
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Figure 5: Spectral function for the one dimensional Hubbard Model within the large U approximation. Upper Panel: quarter- 
filling n = 0.5, computed for U/t = 7.5. Lower Panel: n = 0.7, computed for U/t = 11.5 



nontrivial fermionic matrix elements are computed between wave functions of free fermionic states on a ring, with 
different twisted boundary conditions imposed by the spin configurations. This simpHfication permitted to obtain the 
spectral function in the infinite U limit. Note however that if J7 — > oo the spin spectrum collapses and the spin sector 
is completely degenerate. 

Once the t/U is finite, the problem can be treated perturbatively, and to get the first order corrections of the 
energy it is sufficient to look at the expectation value of the perturbing Hamiltonian (cx t/U) with the unperturbed, 
spin-charge factorized wave functions. When calculating the spectral functions, additional corrections appear in the 
matrix elements that come from applying the unitary transformation to the electron creation and anihilation operators 
[5511^ . For our purposes the most important effect of the finite t/U is to introduce a finite spinon velocity, and that 
is already captured by the first order corrections to the energy. The spinon velocity at the Fermi momenta is given by 



1 



sin 27rn \ 



27112 



O 



1 

IP 



(13) 



where n is the band-filling, and the exponents are calculated at the Fermi level. 

The results of [23 and its extension to finite U were proven to be quite accurate for U/t > 6 (see [50] )• Using 
this method the ID spectral function was obtained considering systems with size L, ranging typically from 120 to 
300 sites; quantitative differences as a function of L were observed to be small within this range. Moreover, in order 
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to reduce the computational time, the results presented here used only contributions from one and two particle-hole 
excitations that were shown to carry the vast majority of the spectral weight (> 99%) jl7]; the inclusion of higher 
order processes was observed to give neglectable contributions. The values of U were obtained fitting the spin velocity 
Vs with the expression (13), after having computed the ID spectral function with an effective exchange constant J^s 
of the order of ~ 0.2. Using the RPA expression ([T]), the 2D spectral function was computed for different values of 
the band filling and transverse momentum. The exact position of the bound state dispersion was obtained as well as 
the new FS and the dependence of the QP weight. The results are presented in the next sections. 



A. Finite large U 

Fig. [sjshows the Hubbard chain [t' = 0) spectral function for quarter-filling [n — 0.5) and a large value of U jt — 7.5 
and for n = 0.7; U/t = 11.5. Close to zero energy (chemical potential) there is a large spectral weight along both the 
spinon and holon branch lines. Note that the spectral weight along the spinon branch dies out as we move away from 
the Fermi level towards positive energies, while the spectral weight along the holon branch line remains high. The 
branch lines for arbitrary values of the Hubbard coupling, [/, are obtained moving one excitation (spinon or holon) 
along its band while keeping the other one fixed at the Fermi level. In the vicinity of the branch line the spectral 
weight has a power law behavior with exponents that may be negative (yielding a large spectral weight) or positive 
(leading to an edge and small spectral weight). As shown in Fig. 1 of Ref. [13J the exponent along the spinon branch 
line for positive energies changes sign from negative to positive and, therefore, there is a loss of spectral weight, while 
the exponent along the holon line is always negative. 

The ID results of Fig. [5] are to be compared with those of Fig. [6] where the spectral function is computed for 
an anisotropic square lattice {t'2 = ^3 = 0), for different values of the interchain coupling and transverse momentum 
and for band-fillings n = 0.5,0.7. Note that since, in this approximation, t'(k) — for ky = 7r/2, the SF for this 
value of the transverse momentum is given by the ID case of Fig. [5] The low energy behavior near the ID Fermi 
momentum agrees with the results of the previous section. In the left panels the effective hopping t'(fc) > and in 
the right panels the effective hopping is t'(fc) < 0. As shown in the previous section this implies that for fcj, < 7r/2 the 
FS increases in size and for ky > it/ 2 the FS shrinks. In the ID case there is a high spectral weight along both the 
spinon and holon branches at the Fermi surface. Introducing the transverse hopping we find as for the coupled LL 
that there is an increased spectral weight in one of the two branches depending on the sign of <'(k): for t'{k) > at 
positive energies the weight is concentrated in the spinon branch and at negative energies in the holon branch while 
the opposite occurs for t'{k) < 0. 

Bound states arise near the spinon branch and their weight increases with t'. For the low energy region the spectral 
weight outside the ID continuum is strictly zero due to phase space constraints. In this case the sharp coherent 
features are poles of the 2D Green's function. Besides the bound states near w = 0, anti-bound states are formed 
at high energies. However for the high energy part of the ID spectrum there is generically no region with strictly 
zero spectral weight since small contributions will come from higher order particle-hole processes not considered in 
our method. This means that in practice, contrarily to the boundstates, anti-boundstates will have a small width 
corresponding to a long, but finite, lifetime of this QP-like excitations. As the transverse hopping increases, all the 
coherent features become sharper inside and outside the ID continuum. However, there is still a significant distribution 
of spectral weight through a continuum, a ID characteristic. Note that a bound state emerges from the edge of the 
Brillouin zone that extends to lower energies, as the transverse hopping grows. 

In Figure [7] we show the 2D Fermi surface and the quasiparticle residues for different transverse hoppings. The 
left panels of Fig. |7] show the evolution of the FS as the interchain hopping is increased. Comparison with the 
non-interacting case (Orange line) shows that interactions will tend to prevent warping of the FS. The Right Panels 
display the value of Z (black lines) and i'(k) (orange lines) along the FS. The QP weight clearly increases with t'^. 
Along the FS the inhomogeneities of Z are quite smooth except for the vicinity of ky = tt/2 where it vanishes. As 
discussed in sec. II A[ since i'(k) = at this point the RPA expression is known to fail. Higher order corrections will 



give a finite Z value leading to a non-zero QP weight along the FS and thus to FL-like behavior. Note also that even 
for <'(k) 7^ the RPA underestimates the value of Z, so higher order corrections will be expected to slightly increase 
its value. 



B. Infinite U 



At infinite U the spinous are dispersionless {JcS = 0) and the spin velocity vanishes — 0. As a consequence, 
the lower edges of the continuum, defined by the spinon dispersion relation, become flat and the continuum in these 
regions extends to zero energy. This is shown in Fig. (Is]). The central panel for ky — tt/2 displays, as before, the 
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ky = ky = tt/4 ky = 3tt/4 ky = tt 




Figure 6: Top 3 rows: Spectral function of the Hubbard model at quarter-filling (n = 0.5) and U = 7.5 in an anisotropic 
square-lattice obtained by weakly coupling Hubbard chains within the RPA ([T]) for different values of the inter-chain hopping 
t{ and transverse momentum ky. The axes labels and the scale are the same as for Figs. [s]. As t'l increases the bound states, 
corresponding to a coherent excitation, start on the boundaries of the continuous region changing the shape of the FS. Lower 
row: Spectral function of the Hubbard model for n = 0.7 and t'l — 0.25 and U/t — 11.5. 

Hubbard chain spectral function. Bound states can still form in the regions where the ID spectral weight is strictly 
zero once t' is introduced. However, this region is smaller than that for the finite U case. Coherent features also appear 
at low energies when the bound states enter the continuum. Considering different values of the transverse momentum 
we see the same trends as for finite U. In the left panels there is a "refraction" of the accumulation of spectral weight 
from a "spinon" branch line at positive energies (note that it is now a fiat line since the spinon velocity vanishes in the 
U ^ oo limit) and a holon branch at negative energies. In the right panels it is the opposite. However, the antibound 
states associated with the holon branch also sharpen, even though the distribution of the spectral weight through the 
continuum is much more visible, as compared to finite U. Since the bound states, associated with the spinous do not 
concentrate much spectral weight, this is to be expected. 
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Figure 7: Left Panel: The Gray and White regions correspond, respectively, to signRe [G (a; = 0, fc^;)] < and 
signRe [G (tj = 0, fco;)] > , i.e. to the exterior and interior of the FS obtained within the RPA. The FS of non-interacting 
fermions is given by the Orange curves. Right Panel: Quasiparticle residue along the FS as function of ky. For each value of 
were determined on each side of the FS. For these values, the bound-state equation was solved in order to find 
UJ+/- — 0; ^+/- (Black lines) was computed using Eq. (Ml. The RPA t' (k) along the FS is plotted as a function of ky (Orange 
curve). As implied by the RPA expression, when the self energy vanishes the coherent excitations disappear [Z = 0). Upper, 
Central and Lower panels correspond respectively to t'l = 0.1; 0.2; 0.5 and to n — 0.5. 



V. THE ROLE OF FRUSTRATION 



It is interesting to see if frustration, in the sence of addition of diagonal terms to the rung ladder-like hoppings, 
has a similar effect of fractionalization in metallic systems as it does in frustrated magnetic systems. In this section 
we investigate the role of frustration in the finite energy behavior of the system comparing a non-frustrated lattice 
(square) with two frustrated lattices, triangular and fully frustrated. 

Fig. |9] shows the spectral function computed for an anisotropic square {t[ = 0.5,t2 — t'^ ~ 0), triangular 
{t[ = 0.25,^2 = 0.25, = 0) and fully frustrated {t[ = 0.05,^2 = -0.2,4 = 0.25) lattices. As the number of frus- 
trated links increases, one observes that the coherent modes are suppressed, as can clearly be seen in Fig. [9] where 
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Figure 8: Spectral function of the Hubbard model for n = 0.5 and t'l = 0.25 for the U = <x where Us — 0. 
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Figure 9: Spectral function of the Hubbard model for several anisotropic lattices: t'l — 0.5, t2 = ^3 = (top row); t'l = 
0.25, t2 = 0.25, = (central row); t'l = 0.05,^2 = — 0.2,^3 = 0.25 (bottom row). Note the decrease of coherent modes and the 
increase of continuum-like features as frustration increases. 



the incoherent continuum, typical from the ID case, carries much more spectral weight when compared with the 
anisotropic square lattice. The reason for the decrease of the coherent features with the degree of frustration, is easy 
to understand at RPA level. The number of bound and anti-bound states due to t' , the spectral weight and the 
distance of the boundstate from the ID continuum all grow with the magnitude of i'(k). Compared to the square 
lattice, the values of t'(k) for frustrated systems vary much more within the Brillouin zone, i.e. even if the maximal 
value of |i'(k)| is the same for both lattices, stronger oscillations are expected for the frustrated case leading to a 
smaller mean value / (ik|<'(k)|, which unfavors the appearance of bound-states. 

In order to give a quantitative measure of the coherent modes we computed the area of the Brillouin zone occupied 
by the bound and anti-bound states (see Fig. 10 ). Starting from a square lattice with t[ ~ 0.2 we have increased 
the total amplitude of the interchain hopping in three different ways. Table |T] shows the evolution of the area of the 
Brillouin zone covered coherent modes. For a square lattice, with a larger t[, one observes a substantial increase of 
the area occupied by the bound and anti-bound states: When the same increment is introduced along there is a 
small decrease of the area and a substantial decrease is observed if ^3 is further increased. 

At low energies, two dimensional spin systems and electronic systems near half-filling (where they can be well 
described hy t — J like models) , are expected to be rather sensitive to frustration and may develop exotic spin- liquid 
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kf, ky TT 




Figure 10: Coherent Modes computed for the anisotropic square lattice with t'l = 0.2 and U/t = 7.5. Left Panels: Spectral 
function computed for ky = and ky = n. Regions with strictly zero spectral weight are shaded in blue, the boundstate (Blue 
line) corresponds to a pole of the 2D Green's function. The high energy regions, shaded in Red, have low but non-vanishing 
spectral weight, therefore the anti-bound states (Red line) arising in this region have a small but finite width. Right Panel: 
Dispersion relation of bound and anti-bound states. 





t'l = 0.2 
= 0.0 
t's = 0.0 




t[ = 0.3 
4 = 0.0 
t's = 0.0 
(+11.0%) 


t[ = 0.2 

t's = 0.0 
(-0.07%) 


t[ = 0.2 
t'2 = 0.1 
t's = 0.1 
(-16.0%) 



Table I: Evolution of the area of the Brillouin zone covered by coherent modes for different values of interchain hopping. The 
percentage values are relative to the area of the square lattice with t'l — 0.2, t'2 = t'^ = 0.0 . 

phases with non-FL behavior. Even if we do not study this low energy regimes, the results presented here do point 
out that the finite energy spectrum is significantly affected by the frustrated nature of the lattice even if the interchain 
hopping is small compared with the monitored energy scale. 

VI. DISCUSSION 

The unusual non-Fermi liquid like properties of some two-dimensional strongly correlated systems has lead to the 
proposal that some signatures of the exotic properties of one-dimensional systems may be observed in their two- 
dimensional counterparts. The dimensional crossover from one to two dimensions has been considered by several 
authors and in most cases it has been found that the one-dimensional features are to a large degree lost, particularly 
at low energies. One characteristic of the one-dimensional systems is the fractionalization of degrees of freedom 
which has, however, been shown to persist in some frustrated magnetic systems via the deconfinement of spinous, 
instead of the coherent magnon-like degrees of freedom characteristic of higher dimensional systems. This apparent 
fractionalization has been confirmed recently as shown, for instance, in [22 . 

In this work we have considered Hubbard chains coupled in non-frustrated and frustrated ways and have studied 
the quasiparticle properties via the spectral function. In order to study the crossover from one to two dimensions 
we considered spatially anisotropic systems where the interchain couplings (hoppings) are small compared to the 
intrachain hoppings. The spectral function of the one-dimensional Hubbard model is in general hard to solve but, in 
some limits, it can be obtained exactly/approximatively such as in the infinite/large U limits. This solution was used 
to obtain, in the RPA, the two-dimensional spectral function. In the low energy regime a small interchain hopping 
leads to the formation of a Fermi surface, as shown before by other authors. The appearance of bound states leads 
to a significant concentration of spectral weight, that extends in some cases to finite energies in a way similar to 
the formation of coherent modes, as expected in a Fermi liquid like system. However a significant weight is also 
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observed spread through a continuum characteristic of fractionalization of degrees of freedom. This is particularly 
found when there is frustration in the hoppings. as evidenced by the increase in spectral weight out of the bound 
states as frustration increases. 

It would be interesting to compare these results with experimental results for anisotropic conductors. However, 
to our knowledge, such systems have not been identified. Some systems show anisotropy but they are not weakly 
coupled, such as the BEDT systems |61i. We expect however that with the advent of fermionic cold atoms in optical 
lattices the predictions of this work may be tested and new classes of exotic two-dimensional systems may be found. 
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Appendix A: Spectral function in restricted Hilbert space 

In this section we discuss the extension of a method, introduced in [22' to study anisotropic anti-ferromagnets, 
to the case of electronic systems. This method permits to write down the 2D spectral function as a function of ID 
quantities, in the limit of small inter-chain coupling. The main ingredient is to restrict the Hilbert space to the 
subspace spanned by eigenstates of decoupled chains with few spinon-chargon pairs. In doing so one neglects some 
processes that are expected to carry low spectral weight. Besides their formal final resemblance, the expression for the 
2D spectral function obtained in this way, follows from fundamentally different approximations than the ones leading 
to the RPA result. However, we will show explicitly that some problems arise when dealing with this approach, that 
lead to inconsistencies that prevented us from applying this method. 

From the exact one dimensional solution of the Hubbard model one finds a multitude of excitations that can be 
identified as coming from charge and spin degrees of freedom. However, for practical purposes, single spin-charge 
excitation characterized by their rapidities carry the vast majority (> 95%) of the spectral weight (see [SO]). From 
small to moderate inter-chain coupling, if no phase transition occurs, such excitations are expected to preserve their 
identity furnishing a natural basis for perturbation theory. The physical picture of the perturbed excitations is given 
by one dimensional fractionalized electron (or hole) that hops coherently between neighboring chains. These two facts: 
small inter-chain coupling and low spectral weight of the other types of excitations, allow significant simplification of 
the problem. The former allows an expansion in small inter-chain coupling and the latter justifies the truncation of 
the Hilbert space to two particle states. 



1. Two-particle states 

Let the ground state of the unperturbed system {t'^ = 0) be denoted by |0) — (gjy |0, y), with |0, y) the ground state 
of the Hubbard Hamiltonian for chain y. Its energy is EoL^Ly, with Eq the mean energy per site and and Ly the 
number of sites in the x and y directions, respectively. From the Bethe Ansatz solution, the two particle states are 
labeled by the charge and spin rapidities , by the value of the component of the spin a and by the total 

charge of the state b = ±1, compared to the GS. Such states can be alternatively labeled by their total energy and 
momentum \ei,kx,<7) — a'^ where £; > and are defined by 

Hid \ei,k^,(T,b) = (e; + EoL^) \ei,k^,a,h) (Al) 
T|£,,fc„a,6) =e''=- \euk^,a,h) (A2) 

where T is the operator that translates the system by one lattice site. For sake of clarity a finite system is considered 
at this stage, the thermodynamic limit being taken only in the final results; therefore, £/ and kx are taken within 
a discrete set of values. The two particle states of the 2D system with momentum k = [k^, ky) are defined as the 
Fourier transform in the y direction of the states with only one excited chain: 

|e,,k,a,6) = -^^e^'^yy \ei,kx,(J,h]y) ®y:^y |0,2/') (A3) 
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By construction these states are orthogonal to each other as well as to the unperturbed GS, |0). The projector to the 
subspace spanned by the two-particle states and the t' = GS is denoted P0+2 = |0) (0| + a i (^ii <j\. 

2. Spectral runction 

The 2D spectral function is obtained as the imaginary part of the retarded Green's function S'pcr (a;,k) = 
— i lim^_j.o IniG^(ti; + zry, k) defined as 

G^(w + ir7,q) 

> 

dt (0| cq,.(i)<.(0) + 4,.(0)cq,.(t) 1 6) 



E 

i,k,cr' ,5 



(6| 



(«) 

k,cr,fc 



(61 ct 



q,cr 



(") 



k.a,& + «^ 



(A4) 



with 1 6) the exact GS of the coupled chains. The second equality was obtained using a complete set of eigenstates 

^jTo-b/ approximated by their projection in 



(n) 



k,<T,& / with energy £'^"^ f, = (^^^k'^i,;, + LxLyE^^. Both |o) and 



the considered subspace and the effective Hamiltonian is given by i/eff = IPo+2-ff]Po+2- 
Using first order perturbation theory in the two-particle subspace one finds: 



|6) ^ |o) 



1 



10) 



= 10) + oft'- 



(A5) 



where the last equality follows since H± acting on |0) creates two electron-like excitations in neighboring chains which 
are out of the subspace. Therefore no corrections to the decoupled GS arise in first order in t' within the considered 
subspace. Since HcS does not couple states with different momentum, total spin or charge one can decompose the 
eigenstates as 



(n) 
k.o-.fc 



E V'k,<T,;,(£;) |k,ei,(T, 6) . 



(A6) 



where the summation index runs only over the ID eigen energies. Computing the matrix elements of H^s one finds 
the Schrodinger equation for the amplitudes 



(A7) 



where 



(A8) 
(A9) 



are pure one dimensional quantities and t'{k.) is the Fourier transform of the transverse hopping matrix ([2]). For 
completeness the ID Green's function in this notation is given by 



b=± I 



Ab^a i£l.,kx)Ab^a jShk^) 

oj + b ei + ir] 



Defining B^^^ — ^6,0- (e/, ^2:) V'k.o-.fcfe/) and using Eqs. (A4A8 A9) the approximate 2D Green's function can 



now be written: 



n,b 



^k,(T,6^k,a,6 

^ + bSE^lb + ^V 



(AlO) 
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which coincides with the ID case when t' = 0. Moreover, the particular form of Eq. ( A7) allows the derivation of the 
following identities: 



1 = t'(k)E 



E 



Ab^a {£l,kx)Ab^a- i£i,kx) 



£l 



(All) 
(A12) 



where the first equality is obtained by simple manipulations of Eq. ( A7 1 and the second follows from imposing unit 
norm to the eigenstates. These equalities enable the definition of the complex valued functions 

^ z-bei 



with the properties 



FlUbSE, 



in) ^ 
\i,a.b) 



(A13) 
(A14) 



So that for a test function p{z), analytic in the vicinity of the real line, one has 



dz p{z) 



1 



b) ^k.cr.b^k,, 



(A15) 
(A16) 



where the contour is taken in the domain of analyticity of p{z) and encircles anti-clockwise all eigen energies bSE^^ j^. 



In particular, using p(z) u-l.+iri ' ^^^^ Green's function (AlOl can be written as 

G^(c. + *r;,k) = 



(A17) 



1 



1 



E 



[F^.M^^)V - t'(k)c^-^ + "/ 
1 



[i^k,a.b(w + «??)] t'(k) 



where the contour does not including the uj + irj pole. 

A remark about this method is in order at this point. Note the RPA expression given by Eq. ([ij so the differences 
between the two approaches can be clearly observed. Contrarily to the RPA it is not possible to define a single analytic 
function F gathering both positive and negative energy contributions. This derives from the fact that in the present 
method Eq.(A17l cannot be given as a function of the ID Green's function. Instead, each branch has to be summed 
separately in Eq. (A17) in order to obtain the same result as in Eq. (AlOl, which is itself a consequence of the fact 
that both 6 = ± sectors are uncoupled by the Schrodinger equation. Care must be taken when the bound states cross 
cj = in ( |A17 1; this would correspond to negative excitation energies arising in the Schrodinger equation, signaling 
an instability of the system. Even though it is still possible to give an expression for the Green's function in this 
case, it would not be physically justified to use this result. Since this happens somewhere in the Brillouin zone for 
the Hubbard model it prevented us to use this method to compute the 2D spectral function. For further comparison 
we give the spectral and the Green's function computed with both methods as a function of the ID spectral function: 
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Restricted Subspace 



G«(c. + zr;,k) = Eb=±i 



[Fk,„,b(tj+i7,)]"'-t'(k) 



where 



Xk,,^.b(£)/'r 



-t'(k)Xk 



,(e)]V[t'(k)xL'.„.,(e)]' 



6 P-I' 



(5 (e - Zp) 



where 



RPA 



G^(0. + Z7?,k) = 



[GlD(w+ijJ,fc,)]-l-t'(k) 



where 



5p(e,k) 



(RPA)" , , , 



[l-t'(k)x<-r''(e)]%[t'(k)x<^r'"(s)]^ 



S{e ~ Zp) 



where 



(RPA)'/ N r, r 1 SpiD(i^.k^) 



In the above expressions the thermodynamic Hmit was taken replacing '^ip{ei,kx,cr,b) by / dv Dk_^^a,b {v) p{v), 
where -Dfc^,CT,b (i') is the one dimensional density of states with quantum numbers kxTU^h, and using the definition 
SpiD kx) = X]& D}^^ rj ii (be) Ai,^„ (be, k^) A(, a- {be, k^). Note that when this replacement is done, the Green's function 
acquires a branch cut in the support of Dk^^tT,b (biy) and coherent contributions from the simple poles for both methods. 
Corrections to the effective Hamiltonian method can be included as in Ref. in (35] by considering a larger subspace 
spanned by states containing higher ordered processes along a chain a nd / o r where more than one chain is in an excited 
state. In the generic case the GS will also have corrections (see Eq. (A5l) of higher order in t' . 



Appendix B: Factorized wave function and spectral function 



In the U — > +0O limit of the Hubbard model the doubly occupied sites are forbidden, and the electrons with opposite 
spins cannot jump over each other - the sequence of the spins of the electrons is fixed. As a consequence, the wave 
functions can be written in a factorized form. 



\^p) = \i'L.,Qi{m^\XN'iQJQ)) (Bi) 
where |x) stands for the spin-part of the N electrons that is defined on a fictitious lattice of the N sites, with the 
wave vector Q = 2nK/N {K = 0, . . . , — 1 is an integer) and fq are some other quantum numbers fSS", '57]. The 
describes the N electrons as spinless free fermions with twisted boundary condition imposed by the spins: 



Lkj — 27r/j 



Q 



where the wave vector Q of the spin wave function appears as a phase shift and Ij — 0, 
and energy of the state are given by 



(B2) 

. , L — 1. The total momentum 



P 



j 



E = -2t^coskj 

j 



Strictly speaking, for U ~ +oo all the spin wave function are degenerate in energy, and the U — > +oo limit is taken 
such that in the ground state the \x) coincides with the ground state of the Heisenberg model, with wave vector 

g = 

The electron addition and removal spectral functions are then given as 



(B3) 



Bik,Lj) = Y,D{Q)BQ{k,Lo) 
Q 



(B4) 
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where the AQ{k,uj) and BQ{k,uj) are coming from the charge, and C{Q) and D{Q) from the spin part of the wave 
function and can be evaluated as described in Ref. ^27). We note the absence of the energy scale in the spin part. 

For finite values oiU/t the spin part gets finite dispersion. As it has been noted in [57] (see also Ref. [SI])) in 
the w-resolved C{Q,uj) and D{Q,uj) the weight is to large extent concentrated on the lower edge of the continuum, 
following the dispersion of the one-spinon branch (we note that the number of spinous in the final states is odd, since 
we had added one spinon to the initial spin wave function), so that 

C{Q,Lo) = C{Q)6{^-es-eQ), 
D{Q,uj) = D{Q)5{u:~es+eQ), 

where Eq is the des Cloizeaux-Pearson dispersion [B3J [M] 

eg = |^cff|sin(Q- V2)| , (B5) 

and 

JeS^nll- — — B6 

\ 2nn J U 

is the effective exchange in the A''-site Heisenberg model of the spin part. After the convolution with the AQ{ijj,k) 
and Bq{(jj, k) charge part, the inclusions of the spinon dispersion given above provides a finite spinon dispersion that 
is seen in Fig. [s] it defines the lower edge of the A(w, k) for the k values between the kp and the ikp, and the lower 
edge of the i?(w, k) for —kp < k < kp. 



Appendix C: RPA and next to leading order corrections 

In this section we rederive the RPA results obtained before by many authors and give explicitly the next to 
leading order corrections. However, since ID correlation functions of higher order are needed in order compute these 
corrections they were not included in the computation of the spectral function in the main text. 

The partition function of the model with Grassmanian sources is given by 

= Z\\ ^gct.t'.c+c+.c+ct.C^^^ (CI) 

where Z|| and (•••)|| are respectively the partition function and the expectation value of an operator in absence of 
interchain coupling and iJ = iJ|| + '^l^'ij'^j- The compact notation c^.t'.c ^ J dr J2k a '^L ai''')^' 0^) ™^ 
(^^ .c = j dr ^jj. ^ „{t)c]i.c,{t) was introduced to improve the readability of the expressions and will be used in the 
rest of this section. Inserting a Grassmanian Hubbard-Stratonovich (HS) field -0 to decouple the hopping term and 
performing a subsequent shift in this field in order to let the term within brackets independent from the sources one 
gets 

Z [C^C] = Z\\ j DiP^Di; e-^^ (C2) 

with 



F = -(t/>UC^).t'-i.(V + C)-ln(e^''=+'='^^ 

(C3) 



= 1 for the physical case, but we will nevertheless perform a saddle point approximation in Eq. (C2) which 
can be seen as an expansion around A oo. This method is similar to the one considered in [Ml- When the HS 
variables are bosons such procedure is equivalent to a given mean field decoupling. The saddle point value is defined 

by ^ (i/^^ + C^) .t'-i + (ct) with (...)|| , = /gV-tc+ctv-V / g^tc+ctv.\ ^ Assuming that t' < 1 the saddle 
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point value is ip = when computed at = 0. To quadratic order one obtains 



1 



l + O 



N 



(C4) 



where S'^Fij — S^-S\ii.F with 4* = (t/j,-!/;^) is the second derivative matrix 



/t'-l + (cct)|| -(CC>|| \ 

which is diagonal since the anomalous terms vanish. The Green's function is obtained taking derivatives with respect 
to the sources 



1 

N "C 

-o 



1 

iV2 



where d^t = <5^t + (^5^t 5qi. stands for the total variation and 5 for explicit one. Using the saddle point condition 
6is, F = Q and total variations of this relation one obtains 



^a,a' ~ ^aa' j^Gm' n'^ an'm' a' + O 

where we have defined the bare {t' = 0) propagator and the propagator at RPA level respectively as 



as well as the four point function 



G\\al' 
Gal' 

= t' 



where stands for the connected correlator. 



(cct), 



al' 



(C5) 

(C6) 
(C7) 



lie. 



(C8) 



In standard notation with k — ky , kj^, cr) and = ^ X^cj X^qo-' expression (C5l translates to 

G{k) = G{k) 



N 



G(g)r4 ik,q) 



-o 



1 



where expressions (C7) and (C8l are respectively given by 



G(fc) = Gil (fc||) ^ -t'(k) 
f^{k,q) = i'(q)rii3(/c||,g||), 



(C9) 



(CIO) 
(Cll) 
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with fc|| = (zwn, k|| , (t) and where Gn is the one-dimensional Green's function. The one dimensional quantity 



cfc|| C911 4||4|| 



ID 



■-m = 



G|| (g||)G|| (fc„)G|| (fc|| 



(C12) 



is given as a function of the ID form factors and propagators. Eqs. ( |C9|C12 ) permit to obtain the 2D propagator 
as a function of the ID quantities only. This expression involves higher order correlation functions for the Hubbard 
chain which are not known at this point. 
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